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>  .Abstract 

In  a  computer  simulation  particles  of  various  energies  were 
followed  along  their  orbits  to  determine  the  probability  of 
their  absorption  during  one  encounter  with  an  inner  satellite 
of  Saturn,  giving  special  consideration  to  leapfrog  and  cork¬ 
screw  escape  factors.  The  time  interval  between  encounters 
may  then  be  used  to  yield  a  measure  of  particle  lifetimes  in 
given  satellite  L  shells.  Energy  and  pitch  angle  dependence 
of  satellite  sweeping  efficiency  is  illustrated.  As  expected, 
the  small  gyroradius  electrons  are  completely  absorbed  in  one 
satellite  encounter,  while  the  large  gyroradius  protons  exhibit 
varying  degrees  of  depletion. 


Introduction 


In  Saturn's  inner  magnetosphere  a  combination  of  trapped  particle 
longitudinal  drift  motion  and  satellite  orbital  motion  will  cause  those 
particles  found  in  a  satellite’s  L  shell  to  regularly  encounter  that 
satellite,  which  absorbs  some  fraction  of  the  particles  with  each 
encounter.  A  particle's  probability  of  absorption  during  one  satellite 
encounter  can  be  combined  with  the  time  between  successive  encounters 
to  estimate  particle  lifetime  in  that  satellite  L  shell,  a  helpful 
parameter  for  any  radial  diffusion  considerations.  Previous  calculations 
along  these  lines  were  done  by  Thomsen  [1977;  1979]  and  Mead  &  Hess 
[1973]  among  others  for  the  satellites  of  Jupiter  —  in  cases  of  a 
"snowplow"  effect  where  the  satellite  completely  clears  its  L  shell 
with  one  sweep,  and  a  "random"  effect  where  non- e quat or i ally  mirroring 
particles  have  a  bounce  period  long  enough  or  relative  drift  fast  enough 
that  they  may  escape  absorption  by  drifting  past  the  satellite  during 
an  extremity  of  their  magnetic  latitude  excursion,  also  known  as 
"leapfrogging"  (see  figure  3) • 

Although  not  considered  at  Jupiter,  the  smaller  diameters  of 
Saturn' 8  satellites  and  the  weaker  magnetic  field  cause  an  additional 
escape  mechanism  to  become  appreciable,  viz,  the  "corkscrew"  effect, 
whereby  the  gyroradius  and  pitch  angle  of  a  particle  allow  its  gyro- 
center  to  pass  within  the  effective  satellite  radius  without  the 
particle  necessarily  being  absorbed.  Absorption  is  assumed  to  take 
place  only  if  the  particle  itself  passes  within  the  geometric  satellite 


radius,  which  will  depend  on  the  phase  of  the  particle's  helical 
path  (figure  l) . 


Figure  1:  (a)  Corkscrew  escape,  particle  not  absorbed, 

(b)  Different  phase,  particle  absorbed. 


-  ■  i 

Sample  comparison  of  factors:  2  MeV  protons  at  Mimas. 


Satellite  radius 
Proton  gyroradius 
Effective  radius 
Bounce  period,  T„ 
Drift  in  Tfi/2  " 


180  km 
262  km 
442  km 


for  a  =  60° 
eq 


30.6  sec  for  a  =  60° 
944  km  eq 


1.  Snowplow  profile  — 


$00  **  Vx>  ua  no  0  ICO  lot  )A)  jy.-o  sec 


All  particles  with  gyrocenters  inside 

of  R  ..  . .  are  removed . 

effective 


Leapfrog  profile  —  Thomsen  [1977]  reduces  the 

2x  R  effective 


absorption  by  a  factor  of  = 
or  more  precisely,  a^fx)  > 


drift  in  lj2 


chord  length  at  x 
drift  in  T_/2 

x5 


The  result  for  2  MeV  protons  «  60°)  at 

Mimas  is  — 


Addition  of  corkscrew  effect  —  refer  to  figure  1, 

drawn  for  2  MeV  protons  (a  *  60®)  at  Mimas.  A 

eQ 

computer  tally  of  the  absorption  of  differently 
phased  particles  yields  the  following  profile  for 
the  combined  effect  of  both  leapfrog  and  corkscrew 
factors  — 


The  other  escape  mechanism  which  was  considered  at  Jupiter  arises 
from  the  inclination  of  a  satellite's  orbit  with  respect  to  the  magnetic 
equator,  allowing  ne ar- equatorial  ly  mirroring  particles  to  slip  by 
during  much  of  the  satellite's  orbit.  This  effect  is  minimal  at  Saturn, 
as  its  inner  satellites  always  remain  very  near  the  magnetic  equator 
and  will  continuously  be  able  to  intercept  even  the  near-equatorially 
mirroring  particles. 

A  aero  offset,  zero  tilt  dipole  magnetic  field  model  is  adopted 
for  the  computations  in  this  report  and  the  inner  satellites  are  then 
assumed  to  remain  on  the  magnetic  equator  (zero  inclination  orbits), 
slight  approximations  which  do  not  significantly  affect  the  outcome. 

The  effect  of  satellite  orbit  eccentricity  on  particle  lifetimes 
is  discussed  in  a  later  section. 

The  profiles  presented  in  this  report  show  the  probability  for  a 
particle  to  escape  absorption  on  one  passing  of  the  satellite,  plotted 
against  the  difference  in  radial  distance  between  the  center  of  the 
satellite  and  the  gyrocenter  of  the  particle  (expressed  in  kilometers). 
While  they  qualitatively  resemble  absorption  "signatures"  it  should 
be  stressed  that  these  profiles  indicate  gyrocenter  depletions .  In  the 
dynamical  situation  those  particles  seen  in  a  satellite  signature  by  an 
omnldirect ional  detector  may  be  displaced  radially  anywhere  up  to  one 
gyroradius  from  their  gyrocenters,  a  smearing  effect,  and  in  practice 
the  satellite  has  of  course  interacted  with  the  particles  many  times, 
so  care  sust  be  taken  in  relating  these  profiles  to  the  observed 
particle  flux  meaaurejnpnts . 


Calculational  Details  and  Discussion 


The  expressions,  parameters  and  constants  used  in  this  analysis 
are  taken  directly  from  Thomsen  and  Van  Allen  [19&0] .  In  particular: 

Dipole  moment  of  Saturn  =  M  =  0.20  gauss  R3  ,  zero  offset,  zero  tilt. 

s 

Keplerian  orbital  angular  velocity  of  satellite 


OJk  = 


4. 19 05  *  IQ ' 
o  3/1 


/a  .  O.OIZ5  \  .  .  r 

t  ^ - J  O-  ‘  *  U tuts  or 


Planetary  rotational  angular  velocity  =  =  I.637  x  K>  rad/sec 

corresponding  to  a  rotational  period  of  10  hrs  39.9  min. 

Averaged  longitudinal  particle  drift  rate  (rotating  reference  frame) 

+  for  protons  ( t  a.$t  ^*.rct )  ,  **  for  e.\er.+roH  5 

+  0.45  -  0.  £>Hk75e*p(-fc. 3^*63.,,,, 

or  fro**  Sc-hrl+z  La*T*rot"ii 

f  _  S.  52 Og  -  O .  3  &  I  ctff  —  p .  b39A  -/h  tVainct 

6  6  •  [l.  38 02  -0.  3h>3  (sm^  t  Vi.n«,t)3 


Particle  angular  velocity  relative  to  satellite  =  ^rel  =  /I  -t  co0  -  oj 

Time  Interval  between  successive  encounters  of  the  particle  with 
the  satellite  *  T£  =  2-rj/ |a)r#J| 

Bounce  period  =  T  =  0.8006  \  sl  >/i  ~  H 

u  [£(E*  2 me Z)J  » 

where  H(a  )  =  I.38  -  0. 32 (sin  a  +  sin1/2®  ) 

eq  eq  eq 


Satellite  parameters 

L 

Diameter 

1979  S2/S4 

2.528 

170  km 

Mimas 

3.092 

360 

Enceladus 

3.968 

600 

Tethys 

4.913 

1040 

Dione 

6.292 

1000 

Rhea 

8.787 

1600 

Table  1:  L  values  from  Thomsen  and  Van  Allen  [1980]. 
Diameters  from  Cruikshank  [1978],  except 
1979  S2/S4  from  Van  Allen,  et.  al.  [1980a] . 


Separating  bovine e  motion  from  drift  motion  the  probability  of  a 
particle's  absorption  in  one  bounce  past  a  satellite  is  computed  as  a 
function  of  distance  between  gyrocenter  and  satellite  center.  For  a 
given  energy  and  pitch  angle  a  family  of  particle  orbits,  each  differing 
in  phase,  is  considered.  The  fraction  of  these  intercepting  the  satellite 
body  is  the  probability  of  such  a  particle's  absorption  (figure  2). 

(100  phases  evenly  spread  over  2rr  radians  are  used  in  ny  computation.) 

Figure  2:  Particles  of  different 
gyrophase,  some  fraction  absorbed. 
Fraction  is  a  function  of  distance 
between  gyrocenter  and  satellite 
center. 


This  probability  is  basically  the  corkscrew  factor  which,  as  suggested 
in  the  introduction,  will  modify  the  leapfrog  profile. 
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A  leapfrog  profile  may  be  computed  in  the  same  manner  by  considering 
a  family  of  bounce  motions  past  a  satellite  (figure  3),  giving  the  same 
results  as  Thomsen's  [1977]  analytic  expression  cited  earlier  in  this 
paper.  (This  would  be  the  final  result  in  the  limit  of  zero  gyroradius.) 


Figure  3:  Plane  of  constant  L.  Gyrocenter  bounce 
motions  a3  particles  drift .  Gyrocenters  intercepting 
satellite  are  removed  in  limit  of  zero  gyroradius. 


To  introduce  the  corkscrew  factor  into  the  leapfrog  picture  the 
gyrocenter  paths  which  pass  within  the  satellite's  effective  radius 
are  weighted  by  the  appropriate  corkscrew  f suitor  (function  of  distance 
from  satellite  center).  Thus  the  curves  shown  in  figure  3,  there 
considered  either  to  be  absorbed  or  to  get  psist,  are  now  assigned 


When  calculated  for  several  values  of  AL  ranging  over  +R  _  , 

-  effective 

from  the  nominal  satellite  L  shell  (figure  5),  the  results  may  be 
plotted  as  a  gyroeenter  absorption  profile  (figure  6),  fraction  of 
particles  remaining  vs.  A.L. 


equatorial  plane 


Figure  5:  Gyroeenter  drift 
path  across  absorption  region. 


Table  2: 

Energy  values  examined  (MeV) 


- - >  cLr  'rf  t 

4LI /  r 

1\  \  . 

[ 

j J  J  So .tell'te  L 

V, 

Reff J 

Electrons 

Protons 

0.1 

0.1 

5*0 

1.0 

0.2 

10.0 

10.0 

0.5 

20.0 

20.0 

1.0 

40.0 

30.0 

2.0 

80.0 

This  was  done  for  electrons  and  protons  of  various  energies  (Table  2) 
at  the  six  satellite  positions  given  in  table  1.  The  plots  each  contain 
results  for  9  different  pitch  angles  (equatorial)  and  are  explained  in 

the  example  of  figure  6.  N 

s (  340  km  \ 
efrawn  fo  Scaler 


frac  ficn  of 
p  «Cr  t  i  c  I «  5 
Swtf'V'i  vi  «a 

sa.t  e  Kite 


Figure  6:  2.0  MeV  protons  at  Mimas,  L  =  3*092,  sat.  radius 

pitch  angles  a  =  10,  20,  30,  40,  50,  60,  70,  80,  89 
eq 


!,?♦! 
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Such  a  plot  clearly  shews  the  progression  of  larger  gyroradii 

(hence  effective  satellite  radius )  with  larger  a  .  (The  omission  of 

eq 

=  90°  avoids  undue  complications  in  the  computer  program.  The 
resulting  profile  would  be  indistinguishible  from  that  of  the  a  =  89° 
substitution.) 

Also  apparent  is  a  sudden  trend  toward  deeper  absorption  at  pitch 
angles  greater  than  70°.  Intuitively  one  can  reason  that  this  sets  in 
as  the  vertical  distance  between  corkscrew  "threads"  becomes  comparable 
with  the  satellite  geometric  diameter  and  the  corkscrew  escape  effect 
is  lost.  For  the  Mimas  example  this  is  verified  in  table  3,  where  the 
distance  in  question  is  listed  under  h.  Note  also  that  the  bounce  time 
is  reduced  at  larger  pitch  angles,  but  the  relative  velocity  is 
correspondingly  increased }  no  there  is  no  appreciable  net  effect  on  the 
leapfrog  escape  mechanism  (see  table  3) • 


a 

eq 

h 

tb 

wrel 

u>  .T^L  R  /2 
rel  B  s' 

■a 

1871  km 

45.2  s 

.000280/s 

1174  km 

■9 

1785 

4l.l 

.000294 

1120 

mm 

1645 

37.7 

.000305 

1068 

ho 

1455 

34.8 

.000316 

1021 

50 

1221 

32.5 

.000325 

979 

60 

950 

30.6 

.000333 

944 

70 

650 

29.2 

.000339 

918 

80 

330 

28.4 

.000343 

902 

89 

33 

28.1 

.000344 

896 

Table  3:  2.0  MeV  protons  at  Mimas,  sat.  diameter  =  360  km. 

Vert.  dist.  between  gyrations  =  h,  affects  corkscrew  factor. 
Del.  drift  in  l/2  bounce  .T  L  R  /2,  leapfrog  effect. 


A  somewhat  puzzling  feature  is  that  hump  growing  in  the  middle 
at  larger  pitch  angles .  Somehow  such  particles  have  a  higher  chance 
of  surviving  by  passing  directly  through  the  center  of  the  satellite. 
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A  clue  is  that  the  feature  begins  to  show  at  a  =  40®  and  becomes 

eq 

more  obvious  as  a  increases .  Looking  at  gyroradius  vs .  pitch  angle 
(Table  4),  =  40°  is  where  the  gyroradius  first  exceeds  the  satellite 

geometric  radius  of  180  km. 


a 

eq 

gyroradius 

10° 

53  km 

20 

103 

30 

151 

4o 

194 

50 

232 

60 

262 

70 

284 

80 

298 

89 

302 

Table  4: 

Gyroradius  vs.  a 
2.0  MeV  protons 
at  Mimas . 


For  a  greater  than  or  equal  to  40°  the  gyrocenter  may  pass  directly 
eq 

through  the  center  of  the  satellite  without  losing  the  particles  to 
absorption.  This  is  a  type  of  corkscrew  escape  effect  independent  of 
the  vertical  thread  spacing,  and  thus  is  not  lost  but  actually  enhanced 
at  greater  pitch  angles. 


So  the  qualitative  features  of  these  profiles  are  readily  explicable. 


A  survey  of  the  full  set  shows: 

1.  Low  energy  electrons  are  subject  to  snowplow  absorption. 

2.  High  energy  electrons  exhibit  the  same  features  as  intermediate 
energy  protons  explained  above. 

3.  Low  energy  protons  give  leapfrog  type  profiles,  little  evidence 
of  corkscrew  escape . 

4.  High  energy  protons  are  virtually  unaffected  by  one  satellite 
passage . 

These  observations  are  compatible  with  the  given  particle ' s  gyroradius 
and  relative  drift  velocity. 

Low  energy  electrons,  for  example,  have  gyroradii  much  less  than 
the  satellite  radius  and  have  a  very  small  relative  drift,  thus  their 
gyrations  and  bounces  are  both  too  small  to  avoid  absorption,  and  the 
satellite  leaves  a  snowplow  type  profile. 

High  energy  protons  with  their  leaps  and  corkscrews  on  the  order 
of  thousands  of  kilometers  will  hardly  notice  a  360  km  satellite.  80  MeV 
protons  at  Mimas  have  gyroradii  up  to  1950  km,  drift  4500  to  5000  km  per 
half  bounce,  and  have  several  thousand  kilometer  spacing  of  their  helix 
threads  (pitch  angle  dependent  of  course).  Only  one  or  two  percent  of 
these  particles  are  absorbed  in  one  satellite  passage .  Therefore  it 
would  take  many  satellite  passages  to  give  the  high  energy  proton 
depletion  seen  by  Pioneer  11  instrumentation  in  the  vicinity  of  Mimas’ 
orbit.  Moreover,  the  radial  diffusion  coefficient  of  these  particles 


L 
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must  be  small  enough  that  the  signature  is  not  appreciably  filled  in 
over  the  time  interval  between  satellite  passages.  This  time  interval, 
given  previously  as  T£  =  2  rf  / 1 I  >  12  minutes  for  80  MeV 

protons  at  Mimas.  With  an  exactly  circular  orbit  it  would  take  Mimas 
N  passages  to  sweep  up  half  of  the  80  MeV  protons,  where  0.98^  =  l/2  ; 

N  =  34.3  >  giving  a  measure  of  particle  lifetime  34.3  x  12  minutes  = 

~4l2  minutes. 

However,  the  slight  eccentricity  of  the  satellites'  orbits  cause 
them  to  sweep  larger  regions  than  merely  -Ref^ect  j^efrom  their  nominal 
L  shell.  As  this  sweeping  region  is  broadened  the  lifetimes  of  particles 
will  increase  proportionately.  E.  g.  between  apo-  and  periapsis  Mimas 
will  undergo  a  radial  excursion  of  ~8000  km  and  so  will  sweep  a  region 
8000  km  +  (2  x  wide.  This  is  ~11500  km  for  80  MeV  protons 

(depending  on  a^),  or  ~  3  times  the  area  Mimas  would  sweep  in  a  zero 
eccentricity  orbit.  Therefore  particles  should  survive  roughly  3  times 
longer  on  the  average.  For  smaller  gyroradius  particles  this  factor 
increases,  reaching  ~  22  (8000km/360km)  as  the  gyroradius  vanishes 

(low  energy  electrons). 

Thomsen  and  Van  Allen  [1980]  noted  that  T_  approaches  infinity  for 
electrons  ~  1  MeV,  the  so-called  resonant  energy.  It  is  speculated  that 
electrons  near  this  energy  would  be  able  to  diffuse  across  the  satellite 
orbit  unaffected  during  the  large  T_  intervals.  [Van  Allen,  et.  al., 

£i 

1980b]  Inwardly  diffusing  electrons  of  other  energies  would  be  swept  up. 
More  work  on  this  fonu  of  energy  filtering  is  in  progress. 
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CTR OCENTER  ABSORPTION  PROFILES 

Electrons  —  at  S2,  Mima s,  Enceladus,  Tethys,  Clone,  and  Rhea 
Protons  —  at  S2,  Mimas,  Enceladus,  Tethys,  Dione,  and  Rhea 
Satellite  diameter  is  indicated  to  scale  at  top  of  each  page. 


Detailed  computer  printout  of  these  profile  values  is  available 

on  request  to  the  author,  listing  additional  parameters  for 

each  case,  each  a  . 

eq 

—  gyro radius,  effective  radius,  resolution  of  profile, 
relative  drift  frequency,  drift  in  1/2  bounce, 
vertical  displacement  in  one  gyration,  gyrofrequency, 
and  drift  period  Tg, 
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